Wave Function Mismatches and Coulomb Drag 
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In this paper I study the topological excitations in a pairing state in double layer systems at 
= 1/2 in the presence of disorders. Due to mismatches between single particle wave functions of 
composite Fermions in different layers, the sensitivity of the Chern number of the pairing composite 
Fermion state, with respect to changes of impurities, is infinity. Consequently, Goldstone mode 
in this pairing state is strongly damped at low temperature. This leads to a unique temperature 
dependence of the drag resistance at low temperature. 



Coulomb drag in double layer systems with Landau 
level filling factor i' ~ 1/2 in each layer has attracted 
much attention for the past few years. When the 

distance between layers d is less than the magnetic length 
Lh, numerical result of [Q shows that the ground state 
can be approximated by a {111} state proposed in 
when d 3> Lh, it was argued that the interlayer attrac- 
tive interaction mediated by the Chern Simons field leads 
to an S-wave pairing BCS state of composite Fermions at 
half filled Landau levels ||^. This is a generalization of 
the idea of having an incompressible quantum Hall liquid 
with even-denominator filling factor discussed some years 
ago |p|-p^ . Attempts have been made to identify these in- 
compressible states proposed in |p|-p^ with pairing com- 
posite Fermion states. As pointed out in ||l^], {331} state 
in 1^, and Pfaffian states in [|ll|,^ can be expressed in 
terms of p-wave pairing composite Fermion states with 
total spin (in the layer index space in the present case) 
equal to zero along different directions; while the hollow- 
core state in [|o| can be connected with a singlet d-wave 
pairing. 

These considerations suggest that low energy transport 
in double layer systems a.t v = 1/2, in general, is deter- 
mined by the properties of a pairing composite Fermion 
state. Especially, in the absence of tunneling between 
two layers and at zero temperature, the Coulomb drag 
in a pairing state is nonzero In contrast, the drag 
between two Fermi liquids vanishes as the temperature 
T goes to zero. The in-plane resistance Rn and the in- 
terlayer drag resistance R12 in double layers in a pairing 
state of composite Fermions can be expressed as H 
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For an S-wave pairing state, Oquasi — 5exp(— A/fcT) is 
the conductivity from thermally excited quasiparticles; 
g — kpl X e'^/h and A is the quasiparticle energy gap of 
the pairing state, or the mobility gap when the quasipar- 
ticles are localized. / is the elastic mean free path of the 
composite Fermions in the presence of disorders, much 
longer than the Fermi wave length kp^. a cond is from 
the composite Fermion condensate, as defined in Eq.l2. 
Meanwhile, 
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Here nf^^^^ are the irreducible compressibility of electrons 
and composite Fermions respectively at wave vector g; 
iVc, e are the composite Fermion superfluid density and 
electron number density. Following Eqs.1,2, Rn = R12 
when composite Fermions have an off-diagonal long range 
order {Ns 7^ 0, a cond = 00) and the double layer system 
becomes incompressible, i.e, k'^(O) = 0. 

In the following I consider the topological excitations 
in a pairing state of composite Fermions in a disordered 
double layer system. Due to mismatches between single 
particle wave functions for composite Fermions in differ- 
ent layers, the pairing state wave function in the presence 
of different impurity configurations is topologically dif- 
ferent. Furthermore, topological excitations have a con- 
tinuous energy spectrum. These low lying topological 
excitations lead to a novel coupled double layer liquid, 
characterized by 
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at low temperature limit. Eq. 3 is valid only in the 
strongly coupled limit when acond > cr quasi ■ 

To obtain these results, let me first consider a double 
layer 2DEG where the impurity scattering potential in 
one layer Vi (r) is different from the impurity potential in 
another layer V2(r). The correlation of these two sets of 
impurity potentials is determined by 
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<Vi{v)V2{v) > 



(4) 



Furthmore, I assume short range impurity potential 
s-i < yi(r)V^i(r') >= s-i < yi(r)Fi(r') >= S{r~r')V^, 
s is the area of the sample. When Vi , V2 are completely 
(un)correlated, 6*12 = (1)0. 

Let me introduce the composite Fermion pairing wave 
function in a general form: A(r,R) =< "01 (R + 
r/2)?/;2(R — r/2) > x(l,2). -01,2 are the composite 
Fermion operators in layer 1,2 and x(lj2) is the pair 
wave function in layer index space. It satisfies the gap 
equation generalized to a double layer system with an 
interlayer attractive interaction constant A12 |14|, 
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A(r,R)= J dr'dr"dR'Ai2(r,r") J de[l - 2ni.(e)] A(r', R') 

[Gjf (R + r'72, R' + r'/2)G^4{Il - r" /2, R' - r' /2) 
+Ggf (R + r'72, R' + r72)G^'^(R - r'72, R' - r'/2)] (5) 

Here npi^) is the Fermi distribution function; 
Gjf'^^(r,r') (Gi^'2'4(r,r')) are the retarded(advanced) 
exact Green function in layer 1 and 2 in the ab- 
sence{presence) of pairing potential. Depending on the 
form of Ai2, the order parameter would have differ- 
ent symmetries. In a clean limit, when Ai2(r, r') = 
\s5{v)6{y'), A(r,R) = 5(r)A,(R)x^(l, 2) corre- 
sponding to an S-wave state. When Ai2(r,r') = 
\k-/V,6{v)\/,5{v'), A(r,R) = V,(5(r)Ap,(R)x^(l, 2) 
representing a p-wave state, i — x,y, x^{^: 2) — x'^(2, 1) 
and x'^(lj2) = — x'^(2,l). Without losing generality, 
I will restrict myself to the S-wave pairing composite 
Fermion state. The final conclusion can be extended to 
a p-wave pairing state. 

Consider the exact Green function at e = and 6a. — 
in 2 D Go^''^( r, r') = cos(fcF|r - r'| - Tr/2 + 0i,2(r)) x 
w/\/fcF|r — r'l, with (/)i,2(r) as the random phase shift at 
point r in the presence of impurity potential Vi,2(j)- By 
the same token, in a given sample, Gj^(r, r')Gg^(r, r') cx 
|r — r'|~^ cos((/>i(r) — (j)2{r)), with the fast oscillating 
phase factor kplr — r'| cancelled out. When Vi{r) is 
not completely correlated with V2(r), (j)i{r) — (/)2(r) is fi- 
nite. GJo^Gqo consists a series of nodes at points r where 
01 (r) — (/)2(r) — {n + l/2)7r, due to mismatches of two 
sets of single particle wave functions at the Fermi energy 
in a double layer system. However, the presence of these 
nodes only reveals itself via an exponential decay of Green 
functions with a given sign when the ensemble averaged is 
taken For example, J de < Gjf (r, r')Gg:l^(r, r') >= 
TO(27r|r — r'p)~^ exp(— 0i2|r — r'\/l), <> represents the 
impurity average. Consequently, following Eq.5, at T = 
the self consistent equation doesn't have a homogeneous 
solution with Sa = and As(r) a constant when 6*12 > 
9c = rAo(^ 1); Aq is the order parameter in the absence 
of disorders and t is the elastic scattering time. 

To account random signs of Green functions in the 
presence of mismatches of wave function, one has to 
take into account higher order contributions in term of 
(A;fO~^; this is beyond the conventional approximation 
used in the theory of disordered superconductors Q. 
Performing an expansion of the self consistent equation 
around 6c diagrammatically (See Fig.l), I obtain for A^, 

leoV' + ^]A.(r) + J rfrif,(r,r')A,(r') = ^A,(r), 

< Ksir,r') >=0, 

< if,(ri,r'i)i^,(r2,r^) > = 

^[Sir, - r,)S{r', - r^)[<5(r, - r',) + J^^^]- (6) 
Here V = V - iSa, 6612 = di2 - Oc, Co = /d/A^. 



Meanwhile, 5a = a -I- A, a is the Chern Simons vector 
potential and A is the vector potential of the external 
magnetic field; W x 6a ~ jA^p — |As„p, which shows 
that the variation of pairing wave function results in 6a. 
(As„ is the uniform solution of Eq.6.) When Kg = 0, 
Eq.6 yields a vortex solution with the size y/6c — 612^0- 
The energy cost of a vortex excitation is gAo{9c — ^12) 
for the paired composite Fermion state coupled with a 
Chern-Simons field. Coulomb energies associated with 
these topological excitations are negligible because they 
extend over a distance much larger than the magnetic 
length. 
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Fig.l. Solid lines are electron Green functions and dashed 
lines are impurity scatterings. 1,2 are layer indices. Di- 
agrams in b) stand for a Hikami Box. Diagrams in c) 
represent the Landau Ginsburg expansion while the dia- 
gram in d) is the typical contribution of Ks{r, r'). 

Eq.6 is a stochastic Landau Ginsburge equation with 
Kg counting spatial random oscillations of G^qCq^ at a 
length scale much larger than l/9i2- factor in the 

variance of Ks{r,r') in Eq.6 reflects the Wigner-Dyson 
level statistics; namely, the fluctuation of number of lev- 
els within an energy band of Thouless energy of size of 
coherence length is unity, regardless of the corresponding 
average number of levels g fl^ . 

The statistical property of solutions of this equation 
was discussed previously in the context of a dirty su- 
perconducting film in Clogston limit p^ , p^ . At 6*12 — 
9c > Sc/g, pairing takes place at a length scale = 
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(,o\/\(^i2 — dc\/&c- The number density of these meso- 
scopic states is l/L'^exp[—g'^\9i2 — 9c\/0c]- Most impor- 
tantly, these mesoscopic states are coupled with exchange 
interactions {J} of random signs. It belongs to the same 
universality class as a XY model with exchange interac- 
tions of random signs. We restrict ourselves to this limit. 

To characterize the ground state and metal stable 
states, I introduce the Chern number of the paired state 
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The constant in front of the integrals is to ensure Mt 
defined above to be equivalent to the standard defini- 
tion of Chern number for the many body wave functions 

>]• 



> 



< 



9* I a* 



A^T - Ady,[< / S^lQ^ 

{xi,yi) are the coordinates of the ith electron. The last 
equality in Eq. 7, which follows Stoke's theory, shows 
that Nt is determined by the phase of the condensate 
wave function along C, the boundary of the sample S. 




Fig. 2. A typical topological excitation. 0, tt, 2tt stand for 
the lift of the condensate phase in those regions. 

First let me point out that Nt of the ground state is a 
nonzero sample specific quantity. Consider a condensate 
wave function, which is topologically different from the 
ground state with A/"^ 7^ Nt- Physically, this new con- 
densate wave function can be obtained by adding a phase 
TT to the condensate wave function in the region enclosed 
by contour Ci , 2tt phase factor to the region enclosed by 
contour C2 while keeping the phase in the region enclosed 
by C3 unchanged. The remaining region is enclosed by 
a contour C4, which has three extremely narrow stripes, 
si,2.3- N^ is the summations of four contour integrals 
defined in Eq.7 along Ci^2,3,4- 

Obviously, the energy cost of the new configuration is 
from singular narrow strips §1^2, along which the phases 
jump by ±7r. Energy of such a configuration is | J|i^/^, 
with i^/^ originating from the summation of random ex- 
change interaction {J} along the stripes. \J\ is a ran- 
dom positive quantity. Let the probability distribution 
function be P(|J|), which is continuous at J = 0. For 



a given L, LP{Jmm)Jrr 



1. Therefore J„ 



is proportional to L}-l'^~^ . 7 > 1/2 means excitations are 
of finite energies, as suggested in JTot . 

Now let me introduce { J' = J + (5J}. The energy cost 
of a configuration with N^ considered above acquires an 
additional contribution of a random sign from the con- 
tour C4, of order of 5JL^^^ . For a given change of ^ J, we 
can always choose a L such that Jmin{L) cx L^'^ <^ SJ. 
In this case, the new ground state belongs to a homotopy 
class different from that of Nt. Excitations of different 
Nt with minima energies of a given homotopy class cor- 
respond to topologically distinct metal stable states. The 
existence of these low energy metal stable states suggests 
possible hysteresis. 

Let e = 1- < Vi(r)V{{r) + FaWFslr') > / < 
Vi{r) + V2^(r) > stand for the correlation between two 
different double layer impurity potentials (Vi,V2) and 
{V(,V^). (< >=< (Vi,2)' >.) Thus, Nt is a 

nonanalytical function of 6, 



lim < 



[Nt{VuV2)-Nt{VI,V{)Y^ 



>= 00. 



(8) 



By contrast, the sensitivity of the Chern number defined 
in Eq.8, is zero for a normal mesoscopic sample, a BCS 
superconductor or a quantum Hall state. 

Infinity sensitivity implies that there are gapless topo- 
logical excitations. To account this, we introduce 



dNr 
dt 
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with as the thermal activation rate of all the possible 
topologically distinct excitations. The thermal activation 
rate for a given barrier E is cT exp{— E / kT) , c is a num- 
ber of order unity; the distribution of the finite energy 
barrier E is assumed to be Pb{E). Naturally, in the pres- 
ence of states nearly degenerating with the ground state, 
Pb(0) > 0. As a resuh, 

-=cT f dVeM-4^)PB{E) « cT^Pb{Q), (10) 

which is a power law function of T instead of an expo- 
nential function at low temperature limit. As far as is 
longer than any microscopic time scale, Eq.9 represents 
the low energy dynamics. 

Consequently, collective modes become strongly 
damped because of thermally excited topologically diff'er- 
ent excitations. The situation here is similar to the gauge 
theory in the hydrodynamic limit proposed for Heisen- 
berg spin glass |Q. The density-density fluctuation and 
current -current fluctuation can be calculated as lEll 



T^dd 



me 



m iLu — Tk 



with 7=1. The typical energy of topological excitations 



which agree with 'Kdd, t^cc obtained in ||^ when topologi- 
cal excitations are absent. Ns and Vs are the superfluid 
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density and the zero sound velocity respectively in the 
absence of topological excitations. The conductivity of 
this state is given as 



0-cond = lim — 7r<id(jw,q) = — 2-w^Tk(T). (12) 

Substituting the results in Eqs.lO, 12 into Eq.l, we ob- 
tain Eq.3. Strictly speaking, at zero temperature 
is nonzero, determined by the quantum tunneling; it is 
however exponentially small when g ^ 1. 

The change of the barrier 5E in the presence of a 
transport current / can be estimated according to min- 
imal coupling principal. 5E = Ti/e J^drj ■ Vx Ih/e, 
with j the current density and Vx the phase gradient 
of the condensate. When SE is of the order of kT, the 
typical energy of topological excitations, the nonlinear 
effect becomes important. This yields a characteristic 
current /„o„ cx kTh/e'^, at which acond is reduced by 
a factor of 1/2 of that given in Eq.l2. It implies that 
i?i2(/ > Inon) < Ruil = 0). On the other hand, at 
this current, the raise of qausiparticle temperature due 
to Jour heat, ST, can be estimated as (Jquasi'Ei'^Tph/Cv, 
with E as the applied electric field, Tph the inelastic scat- 
tering time, Cv the electronic specific heat. It is of order 
of T/g^ X Lpf^/L'^, vanishing in the limit when Lph <C L. 
Here Lph is the inelastic mean free path, and L is the 
sample size. 

Modulations of local chemical potentials due to ran- 
dom impurity potentials lead to pinning of vortices as 
well. For short range random impurity potentials with 
the correlation length Vc ^ ^Oi the pinning potential of 
a vortex is / drSpVi{r) cx gAVo/ep x rc/^o', while the 
vortex excitation energy in a paired composite Fcrmion 
state is gA. When Vo/ep,rc/^Q ^ 1, the pinning ef- 
fect is rather insignificant. For the modulation doped 
GaAs/AlGaAs 2DEG, the correlation length of random 
impurity potentials rc(~ lOOOA) can be much longer than 
the Fermi wave length /c^\- 200 A) and ^o('- 200 A). In 
this case, pinning can be substantial and the coefficient 
c in Eq. 10 should be modified. However, the law in 
Eq.3 remains valid. 

Skyrmions represent a spin texture which is a nontriv- 
ial mapping from a sphere onto itself |^^ . This type ex- 
citation has an energy gap due to Zeeman splitting and 
Coulomb interaction, which is barely affected by disor- 
ders in the present limit. At low temperatures, its con- 
tribution to acond is negligible. 

For a p-wave pairing, Api is suppressed in the pres- 
ence of disorders even when 9i2 — and becomes zero 
when T = Tc ^ I/Aq. Nontheless, one can make the fol- 
lowing substitutions in Eq.6 to obtain an equation for 

T is close to Tc- 5612/ 9c ^ (1 — T/Tc)Sij, 



Api when 



< 



A.- 
with 

Again, Kpij is due to Friedel oscillations with random 



phase shifts, i.e. "mismatches" of wave functions in dif- 
ferent oscillation periods. 

Though a complete microscopic theory for the pairing 
quantum Hall state is still absent, it is conceivable that 
Eq.3 is qualitatively correct for other pairing states such 
as {111} state, {331} state, or Pfaffian states in the pres- 
ence of disorders. Particularly, the stochastic nature of 
Eq.6 is generic to a pairing state involving wave function 
mismatches. In weakly coupled limits the drag conduc- 
tance is related to the topology of the ground state wave 
function, similar to the Hall conductance ||2^]. Recently, 
I. Aleiner showed that the Coulomb drag between two 
disordered Fermi liquids can be of a random sign in a 
mesoscopic limit . 
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